Isotope effect in the presence of a pseudogap 
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We investigate the question whether the unusual doping dependence of the isotope exponent ob- 
served in underdoped high-Tc superconductors might be related to another unusual phenomenon 
observed in these systems: the pseudogap phenomenon. Within different approximations we study 
the influence of a phenomenological pseudogap on the isotope exponent and find that it generally 
strongly increases the isotope exponent, in qualitative agreement with experiments on underdoped 
high- Tc compounds. This result is stable against strong-coupling self-energy corrections and also 
holds for recently proposed spinfluctuation exchange models, if a weak additional electron-phonon 
coupling is considered. 



PACS: 74.72.-h, 74.62.-c, 74.20.-z, 74.20.Mn 
I. INTRODUCTION 

The isotope effect in high- Tc cuprate superconductors 
is unconventional in different respects. Optimally doped 
samples show a very small isotope exponent a of the 
order of 0.05 or even smaller, in contrast to the con- 
ventional Bardeen-Cooper-Schrieffcr (BCS) value of 0.5, 
which one expects .-for a conventional phonon induced 
pairing interaction.EJ This unusually small value in con- 
nection with the high value of Tc lead to early suggestions 
that the pairing interaction in high- Tc cuprates might be 
predominantly electroiiic in origin with a possible small 
phononic contribution.B This scenario, however, is dif- 
ficult to reconcile with the fact that the isotope coeffi- 
cient also shows an unusually strong doping dependence, 
reaching values of 0.5, in some cases |e}ieiji higher, in the 
underdoped, Tc reduced, compounds.BBa 

Many different models have been advanced in order 
to try to understand this unusual doping dependence in 
connection with the small isotope exponent ,at_optimal 
doping, e.g. inftiicnce of van Hove singularitiesErQ, anhar- 
monic phonona3'Ll, electron-phonon coupling iiJ-ttc pres- 
ence of strong autiferromagnetic correlationalBlld, pair 
breaking effectsE.2 —magnetic impurities or Jahn- Teller 
nonadiabaticityt^lEj, but no consensus has been reached 
so far.li3 

In recent years it became apparent that the physics of 
underdoped high- Tc superconductors is governed by the 
pseudogap phenomenon. A behavior which is reminiscent 
of the presence of a pseudogap, growing upon successive 
underdoping, has been observed consistently in a large 
number of different experiments, e.g. nuclear magnetic 
resonance (NMR) Knight-shift and relaxation rate ex- 
periments, specific heat, angular- resolved photoemission 
spectroscopy (ARPES)_iiLnneling, c-axis and a6-plane 
dynamical conductivityllj'Lj. Currently there is no con- 
sensus about the origift|Ot|this pseudogap ^d many dif- 
ferent proposals exist .I13~E3 Williams et al.tZl have shown 
that a phenomenological model for a pseudogap having 
d-wave symmetry can account well for thermodynamic 
quantities in the underdoped cuprates. 

In the present manuscript we want to study whether 



the influence of the pseudogap might give a, perhaps 
more natural, explanation for the unusual isotope effect 
in underdoped iiist-Tc cuprates. As has been shown by 
Carbotte et alHj^cZl an energy dependence of the elec- 
tronic density of states (DOS) varying on the pairing en- 
ergy scale can modify the isotope effect, and therefore 
we ask whether there might be a link between the pseu- 
dogap phenomenon and the isotope effect. Since there 
exists no widely accepted theory for the pseudogap at 
present, here we will follow the idea of Williams et al 
and treat the pseudogap on a phenomenological basis, in- 
troducing it into the single particle excitation spectrum. 
Such a procedure is reasonable, if the pseudogap itself 
does not show aniaotppe effect, as suggested by recent 
NMR experimentsE§N. 

In the following we will study the influence of such a 
pseudogap on the isotope effect within different models. 
We shall start with the weak-coupling (BCS) approxi- 
mation where we consider s- and c?-wave symmetry of 
the superconducting order parameter and of the pseudo- 
gap. In order to see whether these results are stable for 
more realistic cases, we will study two recently proposed 
models based on a spinfluctuation exchange pairing in- 
teraction, e.g. the nearly antiferromagnetic Fenni liquid 
(NAFL) model due to Monthoux and Pinest2l and the 
self-consistent fluctuation-exchange (FLEX) approxima- 
tion for the two dimensional Hubbard model£j 



II. WEAK COUPLING APPROXIMATION 

The linearized gap equation in the weak-coupling limit 
for an anisotropic pairing interaction V{k, k') reads: 



(1) 



Here, is the band dispersion and A(fc) the super- 
conducting gap function. We want to assume that the 
pairing interaction consists of two parts: a phononic 
part Vp{k,k') and an electronic part Ve{k,k'), such that 
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y(fc, k') = Vp{k, k') + Veik, k'). The dominant contribu- 
tion shall be Ve- Then, in weak-coupling approximation 
we have 

V;(fc,fc') = I \^k\,\^k'\ <UJe ^ 

I gIsg 

where tOe is the characteristic energy scale of the elec- 
tronic part and is assumed to be independent of isotopic 
mass, ipvik) is the basis function for the pairing sym- 
metry considered. For s-wave pairing ipsik) — 1, for 
dj.2_y2-wave pairing ^pd ^ 2 (^) — cos20i:/\/2 and for 

dxy-wave pairing ipd^yik) = sin29fc/V2, where 8^ = 
arctan(fcj,/fc2:) is the angular direction of the momentum 
k. 

The phononic part may consist of different contribu- 
tions having different symmetries. However, since we 
assumed that the electronic part is dominating with a 
symmetry specified by ^pr|^ only the VJij-component of Vp, 
having the same symmetry, will affect Tc- Therefore we 
can assume without loss of generality 



Vpik,k') = 



VpQil}n{k)ipri(k') 




if kfel 
else 



(3) 



where ujp is the characteristic phonon energy. In har- 
monic approximation, which we will adopt here, uop varies 
with the isotopic mass M like 1 / \/M, while uje is assumed 
to be independent of M . 

For such an interaction the gap function can be sepa- 
rated into two parts: A(fc) — /S.^{k) + Ap{k), with 



Ae,p(fc) = 



_ J AeO,poV'r,(fc) if kfcl < 







else 



(4) 



With this ansatz Eq. (|^) becomes a 2 x 2 matrix equa- 
tion for the two order parameter components Aeo and 
ApQ. Assuming a cylindrical Fermi surface with a con- 
stant density of states Eq. Im) can be written in the form 



AeO 

Apo 



VeoL{LOe) VeoL{uJe) 
VpoL^We) VpoL{uJp) 



AeO 

Apo 



(5) 



where we defined the function L{uj) 



L{uj) = iV(0) 



tanhfT^^l 
de 



N{0) In 



T 



The last expression holds in the weak coupling limit 
Lu ^ T. N{0) denotes the density of states at the Fermi 
level. In deriving Eq. (||) we assumed < Wp, as is 
usually the case for spinfluctuation exchange models (see 
the following sections). However, the final result Eq. (|l3) 
does not depend on this choice. Letting Lp = L{uip) and 
Lg — L{uje) the leading eigenvalue of the matrix in Eq. 
fH) is 



A(cJe, Wp,T) = —i^ ^ + 



+ ^^(Koie - VpoLp)^ + ^VeoVpoLl (7) 

and Tc is determined from the implicit equation 

X{u!e,UJp,Tc) ^ I. (8) 

From this the isotope exponent ao can be calculated: 



IdlnT^ li 



"0 = n 



2dlnc 



2T 9X dLp d\ PL, 
^ dL„ OTe OLe dT^ 



In the weak coupling limit u>p,uje S> Tc this gives 



"0 = n 



Vpoil-VeoLe 



2 14,0(1 + VeoLe) + Veoil " VpoLp) ' 



(9) 



(10) 



Note, that for a purely electronic interaction Vpo — this 
expression yields ao = and for a purely phononic inter- 
action VeQ = it gives ao — 0.5, as one should expect. 
For a mixed interaction ao will generally lie between 
and 0.5. In fact, one can easily show that for given val- 
ues of ujp and uje one can always choose Vpo and Veo in 
such a way that a given value of Tc and ao G [0, 0.5] is 
reached .e3 

Now we wish to consider the influence of a pseudogap. 
In the presence of a pseudogap we have to modify the 
single particle excitation spectrum. Following Williams 
et al, we replace in Eq. (^ 



(11) 



where Eg{k) is the pseudogap and will be chosen to be 
either Eg s{k) — Eg^ — const for an s-wave pseudogap 
or Eg ,i{k) — Ego cos2Qf^ for a d-wave type pseudogap. 
Note, that this symmetry of the pseudogap not necessar- 
ily has to be identical with the pairing symmetry and we 
will allow them to be independent in this section. How- 
ever, the study in Ref. ^ suggests that both symmetries 
are of d-wave type in underdoped high-T^ compounds and 
we will focus on this case in the following sections. With 
the replacement Eq. ( |ll|) the function L becomes 



2T 



, , 27r uj tanh 
Liio) = ^ / dQ^'i^) / de 



6^ + i?2(e) 



(12) 

Eqs. (^ and (^ still remain valid, if one uses this expres- 
sion for L(uj). In the weak-coupling limit LUp, ujc S> Tc, Eg 
we then find for the isotope exponent 



2t: 



a = ao 



de / de 



^2(9) 



cosh^VfgfM 



(13) 
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where is the isotope exponent Eq. (|lC|) in the absence 
of a pseudogap. Eq. shows that a/ao only depends 
on Ego/Tc, the pairing symmetry V'tvC©) and the sym- 
metry of the pseudogap. Since Tc is a function of EgQ, 
determined from Eq. (j^) , for a given symmetry of both 
the pseudogap and the pairing state a/ao is a universal 
function ofTc/Tco- Here, Tco = Tc{Eg = 0). In Fig. | we 
show a/ao as a function of Tc/Tco for different symme- 
tries. The sohd hne shows the isotope exponent for an s- 
wave pseudogap. This result is independent of the pairing 
symmetry, as can be seen by performing the angular inte- 
gration in Eq. (|l^) . For an anisotropic pseudogap having 
d^2_y2-wa,Ye symmetry, however, the pairing symmetry 
does affect the result. The dotted line shows the result 
for an s-wave superconductor with a dj.2_y2--wave pseu- 
dogap, while the dashed-dotted line shows the result for a 
d2.2„j,2-wave superconductor with a d^2_y2--wsLve pseudo- 
gap. The weakest Tc/Tco dependence is found for a dxy 
superconductor with a dx2-y2-wa.ve pseudogap (dashed 
line). In all cases one can see from Eq. (O) that a/ao 
diverges for Tc — *■ 0. Thus, in principle arbitrarily high 
values of a can be reached. As an illustration experimen- 
tal results on_Pr-doped YBCO are shown in this figure as 
solid squaresJj Here, it should be noted that experimental 
results on different compounds can differ somewhat and 
also vary with the dopant used (see Ref. |l]). Certainly 
these differences need further explanation (e.g. see the 
review in Ref. |l^ and cannot be understood solely due 
to the influence of the pseudogap. Here, we only want to 
focus on the influence of a pseudogap alone and investi- 
gate the general tendency and order of magnitude of the 
effect, which is similar in many compounds. 
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FIG. 1. Weak-coupling result Eq. (jl^) for the isotope ex- 
ponent a/ao as a function of T^/Tco in the presence of a 
pseudogap. ao and Tco denote the values in the absence of 
the pseudogap. The solid line shows the result for an s-wave 
pseudogap. For a d^2_y2-wave pseudogap the results for an 
s-wave pairing symmetry (dotted line), a ci^2_y2-wave pair- 
ing symmetry (dashed-dotted line), and a d^y-wave pairing 
symmetry (dashed line) are shown. The solid squares are ex- 
perimental results on Pr-doped YBCO from Ref. [ll 



As an important conclusion we can draw from these 
weak-coupling results that a pseudogap in general leads 
to an increase of the isotope exponent a over its value ao 
in the absence of a pseudogap. The quantitative size of 
this effect depends on the symmetries of the pseudogap 
and the pairing state. However, the qualitative behavior 
is very similar in all cases. The size of ao can become 
small, if a strong electronic coupling constant Vgo and a 
small phononic coupling Vpo is considered. 

III. STRONG COUPLING EFFECTS: NAFL 
MODEL 

Having seen that a pseudogap can lead to an increase of 
the isotope exponent in a weak-coupling superconductor, 
one might wonder whether this effect will survive in more 
realistic models for superconductivity. In order to see, 
how strong-coupling effects affect the results, we want 
to consider a recently proposed spinfluctuation exchange 
model, the nearly antiferromagnetic_Fermi liquid (NAFL) 
model due to Monthoux and PinesEJ. Within this model 
the pairing interaction is provided by exchange of anti- 
ferromagnetic spinfluctuations and the pairing symmetry 
is d.j;2_y2. The (frequency dependent) pairing interaction 
is given by 

V{q, ivm) = g'^xiq, mn) , (14) 

where 5 is a coupling constant, the Bose-Matsubara 
frequencies, and the spinsusceptibility x is given by 

xiq^ii'm) = ^ (15) 

i + e{q^Q? + y^/^s 

Here, Q — (7r,7r) is the antiferromagnetic wavevector, ^ 
is the magnetic correlation length and bJs the characteris- 
tic spinfluctuation frequency. Using this interaction the 
Migdal-Eliashberg equations for strong-coupling super- 
conductors are solved selfconsistently. Here one has to 
solve for the self-energy E 

T.{k,iujn) = ^ X! ^{f^ - f^' ^^^n - 'ii^n')G{k,iu}n) (16) 

k' .n' 

along with Dyson's equation for the Green's function G 
G{k, luj^) ^—^ (17) 

lUJn - Cfc - ^(k.lUJn) 

selfconsistently. Using this solution, Tc is determined 
from the linearized gap-equation 

(j>{k,iuJn) = -— ^ V{k ~ k' .iiUn - ioJn') 
k' ,n' 

x\G{k,ioj,,)\'<t>{k,iuj^) , (18) 

where (j) is the gap-function. For the bandstructure efc a 
tight-binding band with next nearest neighbor hopping 
has been used in Ref. BO^ and we will adopt that here. 
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In order to have a small nonzero isotope exponent at 
optimal doping we consider coupling to an additional 
phonon mode, the 'buckling' mode studied in Refs. ^0 33 
and This mode provides an attraction in the d.j.2_ 
wave channel and its pairing interaction reads 



Vp{q,ivm) 



Vpoi 



2 1x 



2 



(19) 



We do not expect the main results to depend strongly 
on the details of the electron-phonon spectrum, as long 
as its coupling strength is small compared with the spin- 
fluctuation interaction. It is important, however, that 
the electron-phonon interaction has an attractive compo- 
nent in the d^2_y2-wave channel, as has been discussed 
above Eq. (||). For the calculations we choose the pa- 
rameters given by Monthoux and Pines: ^ — 2.3a, XQ — 
44 states/eV, hopping matrix element t — 250 meV, 
where a is the lattice constant.E2l For the spinfluctua- 
tion frequency LOg we choose three different values O.OSi, 
0.06i, and 0.2t, in order to study the cross-over from 
weak-coupling to strong-coupling behavior. For the char- 
acteristic phonon frequency we choose a typical value of 
ujp — 0.2t. Following Ref. |3^, the interaction strengths 
g and Vpo are adjusted such that the transition tem- 
perature becomes Tco = 90K and the isotope exponent 
reaches ao — 0.05. Results for the coupling constants are 
shown in Table |[ Here, the electron-phonon coupling 
constant Xph is defined in the usual way: 



dn 1 



\h — 2 / — 



1. 



n 



^-Im {Vpiq,n + i5)} (20) 



From Table | we can see that higher coupling constants 
g are required for smaller spinfluctuation frequencies. 
Comparatively small values of the electron-phonon cou- 
pling constant are sufficient to yield an isotope exponent 
ao = 0.05. 

In order to study the influence of a pseudogap we intro- 
duce a c?-wave pseudogap Eg{k) = EgQajs2<dk, as sug- 
gested by the analysis of Williams et al.tZl, into the sin- 
gle particle excitation spectrum by replacing in the single 
particle Green's function Eq. ( p^ ) 



efe + ReS 



±\/(efc+RcS)2 + R2(fc). (21) 





g/t 




Eg ^SUpp 1 TcO 


0.03 


5.1 


0.31 


12.3 


0.06 


3.2 


0.15 


6.5 


0.2 


2.4 


0.10 


4.0 



TABLE I. Coupling constants for the NAFL model with 
an additional coupling to the buckling mode Eq. ( |l9|) for 
different values of the spinfluctuation frequency Us- -Eg.supp 
denotes the value of the pseudogap Ego, which completely 
suppresses Tc- 



Here it is necessary to take into account the real part of 
the self-energy S, since the pseudogap opens at the Fermi 
surface, which is renormalized due to the self-energy. Ta- 
ble I also shows the amplitude of the pseudogap, denoted 
by Eg^supp, which completely suppresses Tc to 0. Ex- 
perimentally, the ratio of i2^_,p to Tco is about 6-15, 
depending on the material.ESO The values found here 
indeed turn out to be of this order of magnitude. Note, 
that the renormalization of the pseudogap Eg{k) due to 
the self-energy E is taken into account in this approxima- 
tion. In contrast to the weak-coupling approximation in 
the previous section, the pseudogap as seen in the density 
of states is washed out now and thus is a real 'pseudo'- 
gapEa 

In Fig. 1^ we show a as a function of Tc for three dif- 
ferent values of the spinfluctuation frequency ujs along 
with the weak-coupling result for a d2,2_j,2-wave super- 
conductor with a d^2_y2-wave pseudogap from Fig. ^. 
For higher values of the spinfluctuation frequency a grad- 
ually approaches the weak-coupling result as one should 
expect. For small lUs = 0.03i there are some deviations 
from the weak-coupling limit. However, the results are 
not affected very much. 
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FIG. 2. Isotope coefiicient q as a function of Tc for the 
NAFL model with an additional coupling to the buckling 
phonon mode. The opening of the pseudogap leads to a sup- 
pression of Tc and an increase of a. The solid line shows the 
corresponding weak-coupling result from Fig. |l]. Results are 
shown for difi'erent values of the characteristic spinfluctuation 
frequency: uis = 0.03/: (dashed line), lJs ~ 0.06/ (dotted line), 
and ujs — 0.2t (dashed-dotted line). In each case the cou- 
pling constants have been adjusted such that Tco = 90K and 
ao = 0.05. 



IV. SELF-CONSISTENT FLEX 
APPROXIMATION 

Within the NAFL model the spinfluctuation pairing 
interaction is flxed and does not change with the elec- 
tronic properties. However, it is clear both experimen- 
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tally and theoretically that the pseudogap does affect the 
spinsusceptibility and thus should affect the spinffuctu- 
ation pairing interaction itself. To study such kind of 
effects it is necessary to treat the electronic properties 
and the electronic pairing interaction in a self-consistent 
way. Such a self-consistent treatment is provided by the. 
so-called fluctuation-exchange (FLEX) approximationEil 
for the two dimensional Hubbard model^nd also yields 
a d^i _y-2-wave superconducting state.t3~El The main dif- 
ference with the NAFL model is that the spinsusceptibil- 
ity is calculated from the interacting Green's functions 
within an RPA-type approximation. Then the pairing 
interaction reads: 



V{q,iVm) = -U- — T--. 7 

2 l~UxQ{q,lVm) 



(22) 



where the bubble susceptibility xo is calculated from 
the fully dressed single-particle Green's function G (Eq. 
(p!?!)) self-consistently: 

(q, iVjn) = X! '^^^ + ^' + i^rn)G{k, iuJn) (23) 
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FIG. 3. Isotope coefficient a as a function of Tc/Tco within 
FLEX approximation with an additional coupling to the buck- 
ling phonon mode. Results for different values of the on-site 
Hubbard repulsion U are shown: U — 3.6t (dotted line), 
U — 3.8t (dashed line), U = 4.0t (dashed-dotted line), and 
U = A.lt (dashed-dot-dotted line). The solid line shows the 
corresponding weak-coupling result from Fig. |l|. 



This guarantees that any change in the one particle 
Green's function is reflected in the pairing interaction and 
vice versa. In Fig. || we show a{Tc/Tco) within FLEX 
approximation with an additional coupling to the buck- 
ling mode already considered in the NAFL model above. 
Results are shown for different values of the on-site Hub- 
bard repulsion U along with the weak-coupling result. 
Here, a bandfiUing of n = 0.84 for a simple tight-binding 
band has been assumed. The electron-phonon coupling 
strength again has been adjusted to give ao = 0.05 for 



Ego = 0. For U < 3.6t, a{Tc/Tco) very much follows 
the weak-coupling limit. Only if U reaches values of the 
order of At or higher, deviations from the weak-coupling 
limit become apparent. For higher values of U, a{Tc/Tco) 
becomes flatter and starts to rise only at smaller val- 
ues of Tc- This is a consequence of the influence of the 
pseudogap on the spinfluctuation pairing interaction. In 
contrast to the NAFL model the opening of the pseu- 
dogap not only reduces Tc via the reduction of single 
particle spectral weight at the Fermi level, but it also 
suppresses the lowest frequency spinfl-uptuations. Since 
these are predominantly pair-breakingi23l23, this increases, 
the effective coupling strength of the spinfluctuationsc2l 
as compared with the phonons and thus leads to a reduc- 
tion of a. 



V. CONCLUSIONS 

We studied the influence of a pseudogap on the iso- 
tope exponent for different models having an electronic 
pairing interaction with a subdominant electron-phonon 
interaction. In the weak-coupling limit we found that the 
introduction of a pseudogap leads to a strong increase 
of the isotope exponent above its value in the absence 
of a pseudogap. For Tc ^ the isotope exponent di- 
verges, allowing arbitrarily high values. The symmetries 
of the order parameter and the pseudogap only lead to 
quantitative, but not qualitative changes of these results. 
Strong-coupling effects within the NAFL model do not 
affect the results very much. The size of the pseudogap 
compared with Tc turns out to be of the right order of 
magnitude. Self-consistent treatment of the spinfluctua- 
tion pairing interaction in the presence of the pseudogap 
can lead to stronger deviations from the weak-coupling 
limit. The general tendency that the isotope exponent 
rises upon opening of the pseudogap still remains, how- 
ever. From these results it does not seem unreasonable 
that the pseudogap indeed can have an important influ- 
ence on the isotope effect and might be, at least partially, 
responsible for the increasing isotope exponent in the un- 
derdoped cuprates. 
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